Mon. Not. R. Astron. Soc. 000, [|-?? (2000) Printed 1 February 2008 (MN 1*1^ style file vl.4) 



Studies of multiple stellar systems — IV. The triple-lined 
spectroscopic system Gliese 644 

Tsevi Mazeh^ David W. Latham^, Elad Goldberg^ 
Guillermo Torres^, Robert P. Stefanik^, Todd J. Henry^'^, 
Shay Zucker\ Orly Gnat^ and Eran O. Ofek^ ^ 

^ School of Physics and Astronomy, Tel Aviv University, Tel Aviv, Israel 

^ Harvard- Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA, 02138, USA 

^Present address: Department of Physics and Astronomy, Georgia State University, 1 Park Place Atlanta, GA 30303-3083 
Submitted July 28 2000; Accepted February 15 2001 

ABSTRACT 

We present a radial-velocity study of the triple-lined system Gliese 644 and derive 
spectroscopic elements for the inner and outer orbits with periods of 2.9655 and 627 
days. We also utilize old visual data, as well as modern speckle and adaptive optics 
observations, to derive a new astrometric solution for the outer orbit. These two orbits 
together allow us to derive masses for each of the three components in the system: 
Ma = 0.410 ± 0.028 (6.9%), Mga = 0.336 ± 0.016 (4.7%), and M^h = 0.304 ± 0.014 
(4.7%) Mq. We suggest that the relative inclination of the two orbits is very small. 
Our individual masses and spectroscopic light ratios for the three M stars in the Gliese 
644 system provide three points for the mass-luminosity relation near the bottom of 
the Main Sequence, where the relation is poorly determined. These three points agree 
well with theoretical models for solar metallicity and an age of 5 Gyr. Our radial 
velocities for Gliese 643 and vB 8, two common-proper- motion companions of Gliese 
644, support the interpretation that all five M stars are moving together in a physically 
bound group. We discuss possible scenarios for the formation and evolution of this 
configuration, such as the formation of all five stars in a sequence of fragmentation 
events leading directly to the hierarchical configuration now observed, versus formation 
in a small N cluster with subsequent dynamical evolution into the present hierarchical 
configuration. 

Key words: methods: data analysis - techniques: radial velocities - stars: binaries: 
spectroscopic - stars: binaries: visual - stars: late-type - stars: individual: Gliese 644, 
Gliese 643, vB 8. 



1 INTRODUCTION 

This paper is the fourth in a series on triple-star systems. 
The overall goal of the series is to contribute to our un- 
derstanding of the formation and evolution of multiple-star 
systems. Paper I (Mazeh, Krymolowski & Latham 1993) 
presented an orbital solution for the single-lined spectro- 
scopic triple star G38-13. Paper II (Krymolowski & Mazeh 
1999) developed an analytical second-order approximation 
for the long-term modulation of the orbital elements of triple 
systems. Paper III (Jha et al. 2000) analyzed the triple- 
lined system HD 109648, presenting observational evidence 

* Some of the observations reported here were obtained with the 
Multiple Mirror Telescope, a joint facility of the Smithsonian In- 
stitution and the University of Arizona. 



for such modulations. The present paper is devoted to the 
nearby triple-lined system Gliese 644. 

Gliese 644 (=Wolf 630=HD 152751=HIP 82817; 
Q=16:55:28.76, (5=-08:20:10.8 [J2000], V = 9.02 mag) is 
a nearby system of M dwarfs at a distance of about 6 pc 
(Gliese 1969). The study of its multiplicity began when 
Kuiper (1934) discovered that Gliese 644 is a visual binary, 
later found to have a period of 1.7 years and semi-major 
axis of 0'.'218 (Voute 1946). Joy (1947) noticed large radial- 
velocity variations for Gliese 644, which led him to suggest 
that one of the two visual components is itself a spectro- 
scopic binary with a period of a few days, making Gliese 
644 one of the nearest triple systems. 

Weis (1982) used photographic plates to derive a pho- 
tocentric orbit and confirmed Fleischer's (1957) suggestion 
that the fainter component of the visual binary, B, is more 
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massive than the primary. He concluded that B, which is 
about 0.1 mag fainter than A in the visual, is the short- 
period spectroscopic binary. Weis derived masses of Ma = 
0.28 and A/b =0.56 Mq, suggesting that the system is com- 
posed of three similar M dwarfs. These masses are consistent 
with the dM2.5 spectral type assigned by Henry, Kirkpatrick 
& Simons (1994) to the blended image of Gliese 644. 

In order to learn more about Gliese 644, we started fif- 
teen years ago to monitor the object spectroscopically. This 
was done within a radial-velocity study of a small sample 
of nearby M dwarfs carried out with the facilities at the 
Harvard-Smithsonian Center for Astrophysics (CfA). From 
the beginning of the project we could see two and occasion- 
ally even three peaks in some of the one-dimensional cross- 
correlation functions, which were obtained by correlating the 
spectra of Gliese 644 against our standard observed M-star 
template. Given the visual orbit of the outer binary, the 
triple- lined nature of Gliese 644 makes this system special, 
as it enables us to derive individual masses for each of the 
three components if spectroscopic solutions can be derived 
for both the inner and the outer orbits. This can add impor- 
tant information about the mass-luminosity relation in the 
solar vicinity (cf. Andersen 1991, Soderhjelm 1999) near the 
bottom of the Main Sequence (e.g., Henry et al. 1999). 

The derivation of radial velocities for all three compo- 
nents in the Gliese 644 system presented an unusually diffi- 
cult challenge, because the inner and outer orbits both have 
small radial-velocity amplitudes, resulting from the nearly 
face-on inclination of the outer orbit, i ~ 165°, from the 
relatively long period of the outer orbit, P ~ 1.7 yr, and 
from the relatively small masses of the three stars, M ~ 0.3 
Mq. The lines from all three stars are rarely resolved in 
our spectra; in most cases all three are blended together. 
To solve the problem of extracting velocities for all three 
stars from blended spectra, we had to wait for the develop- 
ment of a more powerful algorithm than the one-dimensional 
cross-correlation techniques that we had been using. Such an 
algorithm is TODCOR (Zucker & Mazeh 1994), which was 
developed originally as a two-dimensional correlation tech- 
nique for extracting velocities for both stars in double-lined 
spectroscopic binaries, even when the two sets of lines were 
not resolved. The extension of TODCOR to three dimen- 
sions for the analysis of triple-lined systems (Zucker, Torres 
& Mazeh 1995) provided us with the tool that we needed to 
analyze the spectra of Gliese 644. 

Even with TODCOR, the extraction of radial velocities 
for all three components of Gliese 644 proved to be diffi- 
cult. To address this difficulty, we developed a new approach 
for solving spectroscopic orbits for systems with composite 
spectra. The new approach searches a model that includes 
the observed spectra together with the orbital parameters of 
the system. The final model is the one that gives the best 
match between the observed spectra and the corresponding 
set of composite spectra predicted by the model. 

The observations were analyzed with this approach by 
two different teams, one at Tel Aviv and the other at the 
CfA, with some differences in the details at various stages 
of the analysis. The two analyses led to similar results, and 
most of the orbital elements agreed within the internal un- 
certainty estimates. However, the velocity amplitudes for the 
outer orbit disagreed by 3.2 and 3.5 times the internal error 
estimates, suggesting the presence of significant systematic 



difference between the two analyses. Therefore, we present 
the results from both analyses and document the procedures 
used in detail. For the final orbital parameters we adopt sim- 
ply the averages of the values derived by the two analyses. 

Visual observations of the wide pair Gliese 644AB ex- 
tend back more than 50 years. Because some recent un- 
published speckle observations were available to us, we de- 
cided to reanalyze the astrometric orbit of the outer binary. 
Comparison between the results of the spectroscopic and as- 
trometric analyses indicates some possible systematic scale 
differences between the two sets of data. This impression 
was supported by the observations and analysis of Gliese 
644 by Segransan et al. (2001), a work published about 
three months after we submitted our paper, while our paper 
was still being reviewed. Their astrometric orbit, which was 
based solely on speckle and newly obtained adaptive optics 
(AO) measurements and did not include any visual observa- 
tions, yielded a somewhat larger semi-major axis. We find 
a similar trend when we use only the modern astrometric 
observations to estimate the semi-major axis. Because of 
the importance of the scale of the astrometric orbit for the 
mass determinations, and because the visual observations 
are more likely to be afflicted by systematic errors, we de- 
cided to rely only on the modern data to set the scale of the 
orbit. 

The details of the visual analysis are reported in sec- 
tion 2. A new approach for solving spectroscopic orbits of 
systems with composite spectra is described in Section 3, 
together with the details of the Tel Aviv and CfA orbital 
solutions. In Section 4 the solutions for the spectroscopic 
and astrometric orbits are combined to solve for the indi- 
vidual masses of the three M stars in the Gliese 644 system. 
Section 5 discusses in detail the visual and IR photometry 
we have in hand for the A and B components of the sys- 
tem, and Section 6 considers the impact of our results on 
the mass-luminosity relation near the bottom of the Main 
Sequence. Section 7 discusses the relationship of Gliese 644 
to Gliese 643 and vB 8, two nearby common-proper-motion 
companions. Our radial-velocity measurements support the 
interpretation that all five stars are moving together in a 
physically bound system. In the final section we summarize 
our results and discuss the implications for two scenarios de- 
scribing the formation and evolution of binary and multiple 
systems. 



2 THE ASTROMETRIC ORBIT FOR THE 
OUTER BINARY 

2.1 The visual and speckle data 

For nearly 50 years Gliese 644AB was the shortest-period 
(1.7 yr) binary resolved by visual means, with an angular 
separation of ~0'.'2. Only when interferometric and speckle 
techniques became available was it possible to resolve sys- 
tems with even smaller separations. Consequently, Gliese 
644AB attracted a great deal of attention, and a large num- 
ber of visual observations have accumulated since its discov- 
ery. 

The history of the astrometric observations of Gliese 
644AB is far from uniform. During the first five seasons af- 
ter its discovery more than 270 visual measurements were 
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Figure 1. Cumulative distribution of the visual observations. 



Table 1. Unpublished speckle observations for Gliese 644AB. 



Date 


Orbital 


e 


P 




(Year) 


Phase 


(deg) 


(arcsec) 


Type 


1984.3570 


0.797 




0.13 


Ids 


1989.2220 


0.631 




0.17 


Ids 


1990.3530 


0.289 


301 


0.198 


2ds 



made, approximately 200 of them by the same observer, J. 
G. Voute, at the Bosscha Observatory in Lembang, Java. 
Voute was also the first to publish an orbit for the system 
(Votite 1946). In the following 45 years some 80 additional 
visual observations were collected by a number of observers. 
The cumulative distribution of all the visual observations is 
shown in Figure |^, where the irregular pattern is obvious. 

In recent years speckle observations of Gliese 644AB 
have also been made at both visual and infrared wave- 
lengths, providing additional measurements of both the sep- 
aration and the relative brightness of the two wide compo- 
nents. Three of the speckle measurements were obtained by 
one of us (Henry 1991) in the infrared (J, H, K), using the 
Steward Observatory 2.3-m telescope on Kitt Peak, and have 
not appeared in the literature. Two are one-dimensional 
scanning observations in the North-South direction, and one 
is a two-dimensional observation. Descriptions of the instru- 
mentation and observing techniques can be found in Mc- 
Carthy (1986) for the one-dimensional speckle observations, 
and McCarthy et al. (1991) for the two-dimensional speckle 
observation. Finally, Segransan et al. (2001) have very re- 
cently reported five adaptive optics measurements that are 
of very high precision. 

A listing of all the astrometric observations, with the ex- 
ception of the three speckle measurements mentioned above 
and the AO measurements, was provided to us by the late 
Charles E. Worley (U.S. Naval Observatory) and was ex- 
tracted from the Washington Double Star Catalog (Worley 
& Douglass 1996). The three unpublished speckle observa- 
tions are listed in Table ^. The type of measurement is in- 
dicated by "2ds" (two-dimensional speckle) or "Ids" (one- 
dimensional speckle). 



With this wealth of information covering nearly 40 or- 
bital cycles (1934-2000) it should in principle be straightfor- 
ward to obtain a high-quality solution for the visual orbit. 
Indeed, a number of solutions have been published since 
Voute's first determination (e.g., Starikova 1980, Heintz 
1984), although the original orbit is the one listed as "defini- 
tive" in the Fourth Catalog of Orbits of Visual Binary Stars 
(Worley & Heintz 1984). The combination of the old and 
the new data can in principle improve the orbital solution. 
This prompted us to take the present study of Gliese 644 as 
an opportunity to rediscuss the visual orbit. 

2.2 The orbital solution 

Using all the available astrometric observations together 
raises two problematic issues. The first has to do with the 
variety of observers, observing conditions, and techniques 
used to collect the astrometric data since 1934, which could 
give rise to systematic errors. The second and more funda- 
mental issue is associated with the triple nature of Gliese 
644, which introduces tidal interactions between the outer 
and inner orbits. This could produce slow changes in the el- 
ements of the outer orbit (as well as the inner orbit) on 
long timescales (e.g. Mazeh & Shaham 1976, 1977). The 
timescale for a complete cycle of such a modulation is of 
the order of 700 years (Mazeh & Shaham 1979), while the 
actual timescale could be shorter. Thus we cannot rule out a 
priori some small changes in the elements of the wide orbit 
during the 66 years of observational coverage. 

In order to address the first of these issues, we begin 
by pointing out a potentially serious shortcoming of the ob- 
servations obtained by the visual technique. These data are 
dominated in number by the measurements made by Voute, 
who observed with a 60-cm refractor. At visual wavelengths 
the diffraction limit of such an instrument is approximately 
0'.'23, which is essentially the same as the angular separa- 
tion displayed by Gliese 644AB throughout its entire orbit. 
In his original publications Voute reports that the measure- 
ments are frequently only "estimates" made from elongated 
or notched images. Under these conditions angular separa- 
tion measurements by the visual technique must be taken 
with extreme caution, since they have often been shown to 
be biased (see, e.g., Douglass & Worley 1992) and could af- 
fect the semimajor axis of the orbit. Nevertheless, the visual 
observations are potentially valuable for the extended time 
coverage they provide. Consequently, we started by consid- 
ering a solution that incorporates only the visual measure- 
ments and none of the modern speckle or AO observations, 
which could have a different scale. Evidence of just this effect 
is discussed below. 

To compute the orbital solutions described in this sec- 
tion we used standard non-linear least-squares techniques 
based on the Levenberg-Marquardt method (Press et al. 
1992). All position angles were precessed to the epoch 2000. 
The solutions were obtained iteratively, solving for an orbit 
and adopting the root mean square residual as the uncer- 
tainty for the next iteration, until convergence. The errors 
determined through this procedure are 0'.'012 and 0'.'013/p, 
which are surprisingly good for the visual technique. Eight 
observations gave unusually large residuals in these prelim- 
inary solutions (3> 3a"), and were rejected. Observations in 
which only the angular separation was discrepant were not 
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Figure 2. Residuals in 9 and p from an astrometric orbit that 
uses only visual data, shown as a function of position angle. 



entirely rejected, because the position angles can still be ac- 
curate (they are measured essentially independently from 
the separations in both the visual and the speckle tech- 
niques). These observations can contain valuable informa- 
tion on the period of the system, and should be included in 
the solution. 

The residuals from a solution using only the visual data 
are shown in Figure |^ Unusual trends are seen in the resid- 
uals of both 9 and p, plotted here as a function of position 
angle in the orbit. The bottom panel shows an obvious pat- 
tern of parallel curves separated by exactly O'.'Ol, which is 
the precision to which visual observers typically read off the 
separations from their filar micrometers. A less obvious ef- 
fect is seen in the 9 residuals displayed in the top panel, 
which seem to change sign on either side of 90° and 270° 
(shown as vertical dotted lines). In addition, when plotted 
as a function of time (not shown here), the residuals in 9 
are systematically negative by about 5° during 1935, which 
was the first of Voute's very intensive observing seasons for 
Gliese 644. 

To further investigate the latter effect we performed a 
sequence of additional orbital solutions in which we removed 
the older measurements one by one and recomputed the or- 
bital elements in each case. These fits are therefore not in- 
dependent, and perhaps a better way would be to compute 
solutions in a window corresponding to a fixed interval of 
time or a fixed number of observations, sliding it along the 
time axis. Unfortunately the peculiar time history for the 
astrometry of Gliese 644, with more than 75% of the data 
obtained in the initial 10% of the time coverage, does not 
allow this (see Figure |l|). As a result of this exercise, we 
detected sharp changes in the orbital elements as soon as 
the observations in 1935 are excluded. Although tidal ef- 
fects as discussed earlier can produce long-term changes in 



the orbital elements, the discontinuities observed are much 
to abrupt to be real. An alternate explanation would be a 
systematic error in the observations for 1935, the vast major- 
ity of which were made by Voute. From the discussion above 
it appears likely that his position angles include a systematic 
error. This is somewhat unexpected for such a careful and 
experienced observer as Voute, particularly since his later 
observations do not show this effect. There are no clues on 
the source of the systematic error, despite the fact that he 
recorded his observations in great detail, and we can only 
speculate as to the cause. We have therefore chosen to ex- 
clude Voiite's first observing season altogether, after which 
the discontinuities in the orbital elements mentioned above 
disappear completely. 

The remaining visual observations, which cover nearly 
50 years, were used to examine the possibility of changes in 
the orbital elements due to tidal effects. For this we divided 
the measurements into two independent groups, although 
once again as seen in Figure |^ the time sampling is not very 
favorable for this sort of test. As a compromise between 
the number of orbital cycles spanned and the number of 
observations in each group, we chose 1955 as the dividing 
point. Independent solutions using the pre-1955 and post- 
1955 data are given in Table ^ 

The orbital elements derived from the two data sets 
are rather similar, with the exception of the eccentricity. 
The change in e amounts to 3.8 times the combined errors, 
an apparently significant effect. Interestingly enough, this 
is one of the parameters expected to vary due to the tidal 
interaction between the outer and inner orbits, and for some 
triples it displays the most prominent modulation (Mazeh 
& Shaham 1979; Paper II). Although very suggestive, we 
hesitate to place much confidence in this result at the present 
time in view of the systematic trends in the visual data 
illustrated in Figure H, which could alter the shape of the 
orbit in subtle ways. Nevertheless, continued observations 
with modern techniques may well confirm this in the future. 
A combined solution using all the visual data together is 
shown in the third column of the table. 

Because of the potential for systematic errors in the an- 
gular separations of the visual data, a completely indepen- 
dent solution was carried out using only the modern speckle 
and AO measurements. The one-dimensional speckle obser- 
vations also contain useful information on the angular sep- 
aration, and can be included in the fit as well. We assigned 
errors for the (two-dimensional) speckle measurements of 
O'.'Ol in angular separation (p) and O'.'Ol/p in position an- 
gle (9). For the one-dimensional speckle measurements we 
adopted an error of 10 percent. We note, however, that the 
uncertainties reported by Segransan et al. (2001) for their 
five AO measurements are substantially smaller than those 
of the speckle measurements, with the result that they carry 
a very large weight in this new fit. Two of those observations, 
in particular, have quoted errors in the angular separation 
of only 0'.'0005. If allowed to carry this enormous weight 
(400 times that of the speckle technique), these AO mea- 
surements would completely dominate the solution, a risky 
situation when little is known about the systematic errors 
they may be affected by. We have preferred to be conserva- 
tive, and have therefore chosen to set the errors in angular 
separation of these two AO observations to the average of 
the other three, which is 0'.'003. 
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Table 2. Astrometric orbital solutions for Gliese 644AB, excluding the 1935 observations by Votite. 



Element 


Pre-1955 


Post-1955 


All visual 


Modern 


Adopted 




(visual only) 


(visual only) 




(no visual) 




P (yr) 


1.71576 ± 0.00058 


1.71634 ± 0.00088 


1.71687 ± 0.00018 


1.71671 ± 0.00091 


1.717267 ± 0.000039 


P (d) 


626.67 ±0.21 


626.89 ± 0.32 


627.087 ± 0.066 


627.03 ± 0.33 


627.232 ±0.014 


a{") 


0.2169 ±0.0013 


0.2143 ±0.0035 


0.2166 ±0.0013 


0.2238 ± 0.0023 


0.2256 ± 0.0011 


e 


0.0337 ± 0.0026 


0.0666 ± 0.0082 


0.0384 ± 0.0026 


0.0448 ± 0.0066 


0.0433 ± 0.0018 


in 


165.0 ± 1.7 


163.2 ± 5.1 


164.9 ± 1.6 


166.9 ±6.1 


163.1 ± 1.6 


} 


124.4 ± 7.7 


83 ± 17 


118.9 ± 7.2 


98 ± 40 


115 6 + 51 


f22000{°) 


169.7 ±6.4 


140 ± 17 


167.2 ±6.3 


150 ± 26 


163.2 ±3.1 


T 


1988.111 ±0.024 


1988.083 ± 0.036 


1988.129 ±0.017 


1988.122 ±0.067 


1988.143 ±0.011 


g,3 1 p2 (arcsec^yr"^ ) 


0.003468 ± 0.000058 


0.00334 ± 0.00021 


0.003446 ± 0.000063 


0.00380 ± 0.00020 


0.003893 ± 0.000058 


Ng 


231 


39 


270 


14 


284 




219 


37 


258 


13 


13 













1 


1 


Span (yr) 


19.1 


26.0 


48.9 


16.0 


66.0 



The new fit based exclusively on modern observations 
is shown in the fourth column of Table ^ Compared to the 
solution that uses all the visual data, uncertainties in all the 
elements are larger due to the smaller number of observa- 
tions and the shorter time coverage. Still, the elements of 
the two fits are consistent within la, with the exception of 
the semimajor axis {2.7a), which is larger in the modern fit 
than indicated by the visual solution. 

Similar indications come from a semi-independent solu- 
tion published for the astrometric orbit based on data from 
the Hipparcos mission. Gliese 644AB was one of the tar- 
gets observed by the satellite (HIP 82917), and although 
the orbital elements are not actually listed in the Hippar- 
cos Catalogue (ESA 1997), a re-analysis of the intermediate 
transit data was performed by Soderhjelm (1999) as part of 
a study of 205 binary systems. Ground-based observations 
were incorporated to some extent in Soderhjelm's analysis, 
so that his solution for Gliese 644AB is not completely in- 
dependent of ours. Nevertheless, the semimajor axis he ob- 
tained, a = 0'.'23, is larger than suggested by the visual data, 
possibly due to his reliance on the very precise satellite mea- 
surements to set the scale. 

As mentioned earlier, many of the angular separations 
measured visually are merely estimates that were made at 
or under the resolution limit of the telescope, and as such 
they are particularly susceptible to personal (i.e., subjec- 
tive) errors. It seems safer, therefore, to rely on the modern 
speckle and adaptive optics measurements to set the scale, 
since these are considerably more objective in nature. On 
the other hand, the position angles of the visual measure- 
ments may still be useful since they are typically more ac- 
curate than the angular separations (see, e.g., Pannunzio et 
al. 1986). The optimal solution for Gliese 644AB, therefore, 
is one that uses the modern data (both and p) mainly for 
scale, and the visual position angles that provide the time 
coverage to strengthen the orbital period. Subtle effects such 
as those illustrated in the top panel of Figure |^ will tend to 
cancel out over many cycles and will not otherwise affect the 
solution significantly because of the much smaller weight of 
the visual data. The result of this combined fit is given in the 
last column of Table ^ and these are the orbital elements 
we adopt for the remainder of this paper. 
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Figure 3. The astrometric orbit and the observations on the 
plane of the sky. The scale is set solely by speckle and AO mea- 
surements (see text). The plus signs represent the visual obser- 
vations, and the triangles and circles are for the speckle and AO 
measurements, respectively. Motion is clockwise. 

A graphical representation of the orbit and the obser- 
vations on the plane of the sky are shown in Figure ^. As 
expected, the visual observations fall mostly inside the orbit, 
since the angular separations are typically underestimated, 
as discussed above. 

Our new orbital solution is in good agreement with the 
old astrometric solutions, except for the scale. The plane 
of the orbit is highly inclined (only 15° from face on), the 
motion is retrograde, and the orbit is nearly circular, with 
a very small but significant eccentricity. When only an as- 
trometric orbit is available, the values of uj, the longitude 
of the periastron, and Q,, the position angle of the ascend- 
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ing node, are both ambiguous by 180°. The values quoted 
in Table ^ eliminate this ambiguity by taking into account 
the spectroscopic orbit presented below. In particular, the 
relative motion of the secondary is away from the Sun at the 
ascending node. 

The inclination angle we derive, i = 163? 1 ± 1?6, is 
slightly larger than the values by Segransan et al. (2001) 
(i = 160?3 ± 0?1) and Soderhjelm (1999) {i = 161°), al- 
though in the latter case no errors are given and it is difficult 
to evaluate the results. On the other hand, the eccentricity 
found by Soderhjelm (1999), which presumably corresponds 
on average to a fairly recent epoch because of the use of the 
Hipparcos measurements, is e — 0.06, which is very close 
to the value we obtain from the post- 1955 visual data. This 
would seem to support the difference we noted earlier, but 
until stronger evidence is found all we can say at the mo- 
ment is that the ground-based and Hipparcos data cannot 
rule out a possible slight increase in the eccentricity of the 
outer orbit of Gliese 644 on a long timescale. 



3 THE SPECTROSCOPIC ORBITS FOR THE 
INNER AND OUTER BINARIES 

We have monitored the spectrum of Gliese 644 with the 
CfA Digital Speedometers (Latham 1985, 1992), using three 
nearly identical instruments. Most of the spectra were ob- 
tained with the 1.5-m Tillinghast Reflector at the F. L. 
Whipple Observatory atop Mt. Hopkins, Arizona (113 ex- 
posures) and the 1.5-m Wyeth Reflector located at the Oak 
Ridge Observatory in the town of Harvard, Massachusetts 
(77 exposures). A few of the spectra were obtained with 
the MMT, also located at the Whipple Observatory (10 ex- 
posures). Photon-counting intensified Reticon detector sys- 
tems were used to record a single echelle order centered near 
5187 A, with a spectral coverage of 45 A and resolution of 
8.3 kms~^. The signal-to-noise ratio ranged from 7 to 40 
per resolution element in the continuum, with 20 being a 
typical value. The first 194 exposures were obtained over 
a 2645-day span between 1984 and 1991, while the final 6 
exposures were obtained 10 years later in 1997. 

To analyze the data we developed a new approach in 
which a candidate model for the system and its two spec- 
troscopic orbits is judged by how well it reproduces the set 
of observed spectra. For the sake of simplicity we choose 
to discuss our approach only for the case of a double-lined 
binary. Extending our approach to the triple-lined case is 
straightforward. 

3.1 The New Approach 

3.1.1 The old procedure 

In previous studies of double-lined systems (e.g., Mazeh et 
al. 1995, Goldberg et al. 2000) we have derived orbital so- 
lutions in two separate steps. In the flrst step we extracted 
individual radial velocities for the two stars from each of the 
observed spectra. Then in the second step we derived the 
orbital parameters, based on the velocities obtained in the 
first step. The velocity extraction was accomplished using 
TODCOR, a two-dimensional correlation technique (Zucker 
& Mazeh 1994) that assumes the observed spectrum is a 



combination of two spectra, with shifts caused by the radial 
velocities of the two components. The algorithm calculates 
the correlation of the observed spectrum against a combi- 
nation of two templates with different shifts. The result is 
a two-dimensional correlation function, whose maximum is 
expected at the shifts corresponding to the actual veloci- 
ties of the two components. To flnd the orbital solutions of 
the binary systems we used ORB (Paper I) — a code that 
searches the orbital parameter space for a minimum of the 
residuals of the stellar radial velocities. 

To derive the radial velocities of the two stars TOD- 
COR requires two template spectra, one for the primary 
and the other for the secondary. To match the actual com- 
ponents of each binary we chose templates from a library 
of synthetic or observed spectra. To optimize the match be- 
tween the templates and the observed composite spectra, we 
ran TODCOR for a variety of different templates from the 
library. To each pair of templates we assigned a measure of 
its flt to the spectra of the binary - the peak correlation 
obtained by TODCOR for that pair of templates, averaged 
over all the observed spectra. We then chose the pair of tem- 
plates that gave the highest average value for the peak. 

To extract the radial velocities for the two stars, the 
brightness ratio of the two templates, a, is also needed. For 
each observed spectrum TODCOR can derive the radial ve- 
locity of the two components either by accepting a prede- 
termined brightness ratio or by finding the best ratio to fit 
that spectrum at the velocities found. Normally we assumed 
the binary system had the same brightness ratio for all the 
observed spectra, and therefore adopted a single value of a 
for all the spectra. To choose the best brightness ratio we 
either averaged the derived a values over the observed spec- 
tra, or ran TODCOR for a grid of predetermined a values, 
searching for the value that gave the highest average value 
for the peak. 

Choosing the best templates and best value for a com- 
prised the major part of the first step of the old procedure. 
Only after the templates and the brightness ratio had been 
established did we proceed to derive the radial velocities, 
which were then used by ORB to find the orbital solution. 

3.1.2 A modified approach 

The separation of the reduction procedure into the two steps 
described above is somewhat artificial and is not necessar- 
ily the optimum way to flnd the best model for the system. 
We describe here a new approach, based on the philosophy 
that the model as a whole, namely the orbital elements to- 
gether with the choice of templates and brightness ratio, 
should be confronted directly with the observed spectra. In 
short, we propose that for any suggested model of the system 
one should construct a set of predicted spectra and compare 
them with the observed spectra. 

This confrontation can be done in a straightforward 
way. Consider a model consisting of a given set of values 
for the orbital parameters together with a choice of tem- 
plates and brightness ratio. First, derive two stellar veloc- 
ities predicted by the orbital parameters for the timing of 
each of the observations. Second, construct a predicted spec- 
trum for each observation by shifting the templates with the 
predicted velocities and combining the two shifted templates 
with the adopted brightness ratio. To compare the predicted 
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spectrum with the observed one, calculate the correlation 
between the two spectra. Finally, the correlations of the dif- 
ferent observed spectra should be combined into a global 
score, which measures the match of all the spectra to the 
model. 

Note that in this new approach the correlation calcu- 
lated for an individual observed spectrum is the correlation 
at the velocities predicted by the orbital solution. This is 
not necessarily the highest possible correlation for that ob- 
served spectrum. There might be a different pair of velocity 
shifts for the two templates that yields a higher correlation. 
However, in the new approach what counts is the correlation 
corresponding to the two velocities predicted by the specific 
orbital solution, calculated for the time of the observation. 

The best model is the one that yields the highest corre- 
lation score. Therefore, the algorithm we propose performs a 
search for the best model in the complete parameter space, 
which includes the orbital parameters and the possible tem- 
plates and brightness ratios. To perform an efficient search 
we suggest an iterative procedure. The solution found in the 
old approach serves as our first guess. This gives us a choice 
of templates and brightness ratio, together with an orbital 
solution. We then iterate, so that every stage of the iteration 
consists of three steps: 

• A search for the best templates and brightness ratio, 
given the velocities predicted by the previous orbital solu- 
tion. 

• A new TODCOR derivation of stellar velocities, given 
the new templates and brightness ratio. 

• The determination of a new orbital solution with ORB, 
by minimizing the of the residuals with regard to the new 
velocities. 

We require that each step yields a higher correlation 
score, and that the iterations converge. 

The extension of the new approach to triple-lined sys- 
tems is straightforward. We need to search for three tem- 
plates and two brightness ratios, a and (3, one corresponding 
to the ratio between the secondary and the primary bright- 
ness, the other to the ratio between the tertiary and the 
primary. 



3.2 Analysis 

The lines of the three stars of Gliese 644 are blended in 
most of our observed spectra. It is therefore especially im- 
portant to have good template spectra, optimized to match 
each of the three components. In particular, it is impor- 
tant that the template spectra have the correct rotational 
broadening when using TODCOR to determine radial veloc- 
ities from blended spectra. If there is not enough rotational 
broadening in the templates, TODCOR is forced to pick or- 
bital velocities that are too far apart, in order to match the 
observed line broadening. Conversely, TODCOR will pick 
orbital velocities that are too close if the templates have too 
much rotational broadening. 

For most of our radial- velocity projects we use template 
spectra drawn from an extensive library of synthetic spectra 
calculated by Jon Morse for a large grid of Kurucz model 
atmospheres (cf. Nordstrom et al. 1994, Morse & Kurucz in 
preparation). Synthetic templates have the advantage that 
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Table 3. The stars observed for templates. 



Gliese a (J2000) 5 V Spectral 

(mag) Type 



OoU 


1 n-1 1 -99 1 A 


-1-4Q-97-1 ^ '\ 




K7 


ouy 


90- ^'^■l Q 7Q 


-i-fi9-nQ- 1 ^ R 
T^uz.uy.-Lo.o 




MO 


846 


22:02:10.27 


-1-01:24:00.8 


9.16 


M0.5 


908 


23:49:12.53 


-(-02:24:04.4 


8.98 


Ml.O 


15A 


00:18:22.89 


-1-44:01:22.6 


8.07 


Ml. 5 


49 


01:02:38.87 


+62:20:42.2 


9.56 


M2.0 


745B 


19:07:13. 


4-20:52:36. 


10.75 


M2.0 


48 


01:02:32.23 


-1-71:40:47.3 


9.96 


M2.5 


725A 


18:42:46.66 


-1-59:37:50.0 


8.91 


M3.0 


725B 


18:42:46.90 


-1-59:37:36.6 


9.69 


M3.5 


273 


07:27:24.50 


+05:13:32.8 


9.89 


M3.5 


699 


17:57:48.50 


+04:41:36.2 


9.54 


M4.0 


51 


01:03:12. 


+62:21:54. 


13.66 


M5.0 


905 


23:41:54.0 


+44:09:32. 


12.28 


M5.5 



they can be calculated for dense grids in effective temper- 
ature, surface gravity, metallicity, and rotational broaden- 
ing. If observed spectra are used as templates, it is hard to 
find real stars that fill densely such grids. Unfortunately, the 
Kurucz models start to become increasingly unrealistic for 
effective temperatures below about 4000 K. For M dwarfs 
we have found that we get better correlations if we use tem- 
plates drawn from a small library of observed spectra with 
spectral types in the range K7 to M5.5. The observed tem- 
plates available for our analysis of Gliese 644 are summarized 
in Table ^ where we list the stellar spectral types of Henry 
et al. (1994). 

Throughout our iterations two sets of templates per- 
formed significantly better than any others we tried: Gliese 
725A for the primary and Gliese 725B for the secondary and 
tertiary, and Gliese 725A for the primary and Gliese 273 
for the secondary and tertiary. These sets gave essentially 
the same peak correlation averages. We decided to adopt 
the Gliese 725 templates, because we hoped that this would 
minimize systematic errors due to differences in velocity zero 
points, metallicity, and rotational velocity. 

To find the brightness ratio and the orbital parameters 
we have developed two independent procedures, one at Tel 
Aviv and the other at the CfA. Completely independent 
codes were used throughout, both for the implementation 
of TODCOR and for the orbital solutions. We describe the 
two procedures in the next two subsections. 



3.2.1 The Tel Aviv Analysts 

To find the brightness ratios of Gliese 644 we ran in each 
iteration a grid of a and (3, using all 200 velocities derived 
in the previous iteration. To average the different correla- 
tions calculated for the different observed spectra we used a 
"generalized correlation score" derived by Zucker (in prepa- 
ration) by applying a maximum-likelihood approach. This 
score can be thought of as a way of weighting the individ- 
ual correlations by using the estimated signal-to-noise ra- 
tio of each spectrum. We finally converged to a value of 
Q = 0.56 ± 0.06 and /3 = 0.36 ± 0.04. 

Finding the orbital solutions of a triple-lined system for 
a given set of velocities is somewhat more complicated than 
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in the case of a double-lined binary. This is so because we 
have to solve for two orbits, one for the close pair and the 
other for the center of mass of the close pair together with 
the distant third star. To solve simultaneously for the two 
orbits we use the code ORB. 20, a slightly modified version 
of ORB. 18 (Paper I). The solution is reached in a few itera- 
tive steps. First, ORB solves for the orbital elements of the 
primary, whose sole motion is within the wide orbit. Second, 
ORB solves for the orbital elements of the primary and the 
secondary as a double-lined system in the wide orbit, ignor- 
ing at this stage the short-period motion of the secondary. 
Third, it solves for the elements of the primary and the sec- 
ondary as a triple system, taking into account for the first 
time the secondary short-period motion. At this stage the 
wide orbit is considered as a double- lined orbit, while the 
short-period orbit is considered as a single-lined one. Then, 
in the last stage of the solution, the tertiary data are taken 
into account, and the short-period motion is considered as 
a double-lined orbit. Each step of the iteration is used as 
the starting point for the next one. Finally, we compute dif- 
ferential corrections for all the elements of the two orbits 
together. 

It has been our experience that the period is the critical 
parameter of the orbital solution, so ORB first searches the 
orbital parameter space with a dense grid of fixed values of 
the orbital period, covering a range of periods given by the 
user. For every value of the period ORB searches for the 
parameters that minimize . This procedure ensures that 
we do not miss any local minimum in the parameter space 
of the orbital elements. Only when the best period has been 
found, does ORB execute an iterative algorithm to find the 
best set of parameters to fit the data. 

The velocities derived with TODCOR can encounter 
two problems, both resulting from the fact that all observed 
spectra include some noise. The problems have to do with 
the fact that TODCOR derives the three velocities from 
the location of the highest peak of the correlation in the 
three-dimensional velocity space. One problem occurs when 
a spurious peak of the three-dimensional correlation gets 
randomly enhanced by the noise and becomes higher than 
the peak that corresponds to the actual three velocities. In 
such a case TODCOR can chose a completely wrong peak. 
Another problem occurs when TODCOR switches between 
the velocities of the three stars. This can happen when the 
templates are similar and the spectra have poor signal-to- 
noise ratios. Both problems can generate outlier velocities 
with large residuals. 

To minimize the problem of identifying a completely 
wrong peak we let TODCOR search for a peak in a lim- 
ited region of the three dimensional space, centered on the 
point which corresponds to the three velocities predicted by 
the orbital solution from the previous iteration. Usually a 
range of ±35 kms~^ in each dimension is searched. To iden- 
tify velocities switched by TODCOR we search for exposures 
where the three derived velocities yield a high when com- 
pared to the velocities predicted by the orbital solution. We 
then switch the velocities if and only if the switched veloc- 
ities give a significantly smaller for that exposure. After 
switching velocities we then iterate the entire orbital solu- 
tion again and require that the new solution give a lower 



The individual radial velocities and (0-C) residuals 
from the orbital solutions are reported in Table 

The elements for the spectroscopic orbits derived by the 
Tel Aviv team are given in Table ^ For four of the observed 
spectra the assignment of the three velocities to the three 
stars in the system was swapped from the initial assignment 
made by TODCOR in order to improve the overall of 
the solution, as described above. We denote the orbit of 
Ba and Bb by B, and the orbit of A and B by AB. The 
period, eccentricity, longitude of periastron, and the time of 
periastron passage are denoted by P, e, u) and T, respectively. 
The radial- velocity semi-amplitudes of A and B are denoted 
by Ka and A'b. Similar notations are used for Ba and Bb. 
We also report the light ratios derived using TODCOR in 
Table |^. The errors in the light ratios were estimated using 
an analysis of the values near the peak of the three- 
dimensional correlation surface. These errors should include 
the effects of photon noise, but do not take into account 
possible systematic errors due to template mismatch. 

Figure ^ depicts the radial- velocity curve for the inner 
orbit of the close spectroscopic pair Ba and Bb according 
to the Tel Aviv solution. The velocities observed for Ba and 
Bb are plotted with the motion of the center of mass of the 
B system around A removed. The radial- velocity curve for 
the outer orbit of the wide visual pair A and B is presented 
in Figure ^. The velocity plotted for B is the average of the 
velocities observed for Ba and Bb, weighted according to 
their mass ratio — i^Ba : i^Bb. 

3.2.2 The CfA analysis 

At CfA the standard approach has been to select the opti- 
mum templates and determine the light ratio for a double- 
lined binary by working with a subset of the observed spec- 
tra, just the ones where the lines of the two stars were well 
resolved. This did not work well for Gliese 644, because the 
lines were resolved in only a few spectra. We addressed this 
problem by sorting the spectra into an order ranked accord- 
ing to the minimum velocity separation between any two of 
the components in the system. Then we calculated the aver- 
age light ratios for all the spectra with velocity separations 
larger than some specified minimum value. Plots of the av- 
erage light ratio versus minimum separation were used to 
select the light ratios for the next iteration. Figure ^ shows 
these a and (3 plots for the final CfA solution. The fact that 
these plots are quite fiat is an indication that we have ar- 
rived at a robust solution where the light ratio is essentially 
independent of blending, down to a minimum separation of 
5 kms~^. For each point plotted in Figure |^ the error bars 
are the standard deviation of the mean value of the light 
ratio. The scatter of the points is similar to the error bars, 
suggesting that these errors are reasonable estimates of the 
internal precision. The errors quoted for the CfA light ratios 
in Table ^ are estimated from these error bars and scatter, 
and are much smaller than the corresponding errors reported 
for the Tel Aviv analysis. 

For each iteration in the CfA analysis we used the or- 
bital solution from the previous iteration to provide TOD- 
COR with initial values for the velocities of the three stars. 
We then allowed TODCOR to search for the nearest peak 
in correlation space. This procedure should eliminate or at 
least minimize the cases where TODCOR finds the wrong 
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Table 4. Radial velocities and residuals (km s~^) for Gliese 644 (first 20 lines) 



HJD Gliese 644A Gliese 644Ba Gliese 644Bb 

(2,400,000+) CfA Tel Aviv CfA Tel Aviv CfA Tel Aviv 

Vr (O - C) Vr (O - C) Vr (O - C) Vr (O - C) Vr (O - C) Vr (O - C) 



45817.8613 


10, 


.88 


— 0. 


42 


11, 


,84 


+0, 


,80 


30, 


.70 


+8, 


,90 


30. 


22 


— 0, 


.32 


11 


.52 


— 0, 


,68 


6 


.00 


+2, 


94 


46461.9627 


14, 


.50 


— 1, 


12 


11, 


,81 


+0, 


.21 


10, 


.80 


+4, 


,02 


16. 


38 


+4, 


,50 


27, 


.24 


+2, 


,46 


26 


.34 


+3. 


50 


46464.9225 


11 


.00 


+3. 


89 


14, 


.55 


+2, 


.84 


15. 


.95 


+2, 


,93 


9. 


33 


—2, 


.71 


25 


.86 


— 1, 


.14 


24, 


.62 


+2, 


11 


46485.9236 


12 


.21 


— 1, 


47 


12, 


.13 


— 0, 


.40 


2 


.44 


— 0, 


.15 


2, 


91 


— 0, 


.93 


30, 


.81 


— 0, 


.69 


29 


.15 


— 1, 


25 


46487.9269 


13 


.30 


— 0. 


21 


13, 


.01 


+0, 


.39 


32, 


.38 


— 1, 


.36 


32, 


42 


+0, 


.66 


—2 


.20 


+0, 


.32 


1 


.08 


+1. 


59 


46489.9227 


13 


.79 


— 1, 


94 


11, 


,69 


— 1, 


.01 


13, 


.02 


+3, 


,03 


16. 


63 


+1, 


,66 


21 


.57 


+0, 


,73 


18 


.04 


+0. 


16 


46492.9182 


12 


.72 


+0. 


18 


11, 


,60 


— 1, 


.24 


16, 


.14 


+1. 


,22 


13. 


09 


—2, 


,83 


18 


.49 


— 0, 


,47 


16 


.24 


— 0. 


40 


46511.9706 


14, 


.13 


+0. 


,45 


13, 


.66 


-0, 


.06 


25, 


,44 


-1, 


.26 


23. 


,49 


-0, 


.87 


4, 


.91 


+0, 


,07 


5, 


.38 


-0. 


,71 


46512.8394 


14, 


.77 


+1. 


,37 


13, 


.92 


+0, 


.16 


0, 


.42 


+0. 


,02 


-0. 


,22 


+0, 


.88 


34, 


.80 


+0, 


,67 


31, 


.27 


-2, 


,85 


46513.8548 


14, 


.82 


-0. 


15 


14, 


,01 


+0, 


,20 


21, 


,46 


-1, 


,52 


21. 


,59 


+0, 


.25 


8, 


.26 


+0, 


,67 


7, 


.63 


-1, 


,66 


46520.0076 


14, 


.34 


-0. 


,71 


14, 


,01 


-0, 


.10 


27, 


.21 


-1, 


,20 


25. 


,54 


-1, 


.79 


1, 


.23 


-0, 


,12 


2, 


.75 


+0. 


,49 


46520.8817 


15, 


.08 


+0. 


,03 


14, 


.93 


+0, 


.78 


24, 


.31 


-0, 


.47 


23. 


,74 


+0, 


.10 


5, 


.93 


+0. 


.58 


6, 


.23 


-0. 


,04 


46523.8260 


14, 


.33 


+0. 


,45 


13, 


,75 


-0, 


.54 


25, 


.28 


-4, 


,27 


22. 


,28 


-1, 


.90 


1, 


.32 


-0, 


,32 


4, 


.19 


-1, 


,29 


46537.9235 


15, 


.75 


+0. 


,79 


15, 


,43 


+0, 


.43 


30, 


,70 


+0. 


,90 


28. 


,82 


-0, 


.29 


-0, 


.04 


+0, 


,38 


0, 


.46 


+1. 


,42 


46538.9682 


12, 


.11 


+5. 


,30 


12, 


,37 


-2, 


,68 


18, 


,73 


+3, 


,45 


18. 


,90 


+5, 


,80 


20, 


.65 


-3, 


,31 


17, 


.43 


+0. 


,80 


46539.8366 


15, 


.35 


-0. 


,43 


15, 


.50 


+0, 


.41 


-1, 


.29 


-0, 


,54 


-1, 


,39 


-0, 


.38 


32, 


.39 


-0, 


.12 


31, 


.86 


~0. 


,28 


46540.8298 


15, 


.00 


-0. 


,40 


15, 


.75 


+0, 


.61 


28, 


.30 


-0, 


.80 


27. 


17 


-0, 


.81 


-0, 


.60 


-0, 


,52 


1, 


.34 


+1. 


,24 


46540.9278 


15, 


.31 


-0. 


,24 


15, 


,22 


+0, 


.07 


30, 


.19 


-0, 


,65 


29. 


13 


-0, 


,46 


-2, 


.37 


-0, 


,21 


-3 


.06 


-1, 


,36 


46541.8861 


14, 


.86 


+1. 


,45 


14, 


,03 


-1, 


.17 


16, 


.54 


-3, 


,58 


13. 


,08 


-1, 


.63 


11, 


.59 


-0, 


,71 


17, 


.12 


+2, 


,48 


46565.7149 


16, 


.63 


+1. 


18 


15, 


.66 


-0, 


.72 


11, 


.68 


+1. 


.31 


11, 


,81 


+1, 


.66 


20, 


.00 


-0, 


.15 


20, 


.92 


+2, 


,91 



Table 5. The spectroscopic elements for the inner and outer orbits. 



Element 


Tel Aviv 


CfA 


Adopted 


Pb (days) 


2.965530 ± 0.000011 


2.965515 ±0.000018 


2.965522 ± 0.000014 


Kb^ (kms-i) 


16.81 ±0.14 


17.20 ±0.17 


17.01 ± 0.20 


K^h (kms-i) 


18.54 ±0.19 


19.00 ±0.21 


18.77 ±0.23 


eb 


0.030 ± 0.006 


0.021 ± 0.007 


0.026 ± 0.007 


^b(°) 


170 ± 12 


162 ±21 


166 ± 16 


Tb (HJD-2440000) 


7337.3 ±0.1 


7337.23 ±0.18 


7337.3 ± 0.14 


OBaSiniB (Gm) 


0.685 ± 0.006 


0.701 ±0.007 


0.693 ± 0.008 


OBbSiniB (GM) 


0.756 ±0.008 


0.775 ± 0.009 


0.766 ±0.009 


MBa Sin3 IB (Mo) 


0.00713 ± 0.00016 


0.00765 ±0.00019 


0.00739 ± 0.00026 


MBb sin3 iB(M0) 


0.00646 ± 0.00013 


0.00692 ± 0.00016 


0.00669 ± 0.00023 


qs = MBb /MBa 


0.906 ±0.012 


0.905 ±0.014 


0.906 ±0.015 


CTBa (kms-1) 


1.35 


1.7 




CTBb (kms-1) 


1.88 


2.1 




Pab (days) 


625.9 ± 1.1 


625.7 ± 1.6 


625.8 ± 1.3 


Ka (kms-l) 


4.81 ±0.08 


4.38 ±0.11 


4.60 ±0.22 


Kb (kms-i) 


3.22 ±0.10 


2.67 ±0.12 


2.95 ±0.28 


fiAB 


0.025 ±0.016 


0.080 ± 0.024 


0.053 ±0.028 


t^AB(°) 


311 ± 36 


285 ± 17 


298 ± 27 


Tab (HJD-2440000) 


7232 ± 62 


7185 ± 28 


7208 ± 45 


OA sin JAB (Gm) 


41.42 ± 0.67 


37.53 ±0.96 


39.5 ±2.0 


aB siniAB (Gm) 


27.72 ± 0.87 


22.93 ± 1.08 


25.3 ±2.4 


MAsin3 jAB(M0) 


0.01348 ± 0.00080 


0.00853 ± 0.00075 


0.0110 ±0.0025 


MBsin^ jAB(M0) 


0.02014 ± 0.00089 


0.0140 ± 0.0009 


0.0171 ±0.0031 


(JAB = Mb /Ma 


1.494 ±0.048 


1.636 ±0.087 


1.56 ±0.17 


CTA (kms"-') 


0.84 


1.2 




7 (kms"-*-) 


14.89 ± 0.05 


15.30 ±0.06 


15.10 ±0.21 


Nab, 


200 X 3 


200 X 3 




a 


0.56 ±0.06 


0.571 ±0.008 


0.566 ±0.034 


/3 


0.36 ± 0.04 


0.355 ± 0.007 


0.358 ±0.024 



© 2000 RAS, MNRAS 000, 



10 T. Mazeh et al. 



.62 

< 6 
m 

II 

B -58 
o 

I .56 
•jj .54 
.52 



J , , , , \ , , , , u 

5 10 16 

Velocity Separation {km/&) 



■n 



J , , , , \ , , , , u 

5 10 15 

Velocity Separation (km/s) 



Figure 6. The light ratios a and /3 from the CfA analysis as a function of the minimum velocity separation between any two components 
of the Gliese 644 system. 
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Figure 4. The radial-velocity curve for the inner orbit. The mo- 
tion of the center of mass of the B system around A has been 
removed from the velocities plotted for Ba and Bb, represented 
by filled triangles and open squares, respectively. 



Figure 5. The radial-velocity curve for the outer orbit. The ve- 
locities observed for A are plotted as filled circles, and open dia- 
monds are used for the average velocity of Ba and Bb, weighted 
according to their mass ratio. 



correlation peak and assigns an incorrect velocity to one or 
more of the stars. Therefore, for the CfA analysis we chose 
not to swap any of TODCOR's velocity assignments. The 
orbital solution was carried out with an independent code 
that solves simultaneously for all elements of the inner and 
outer orbits, by computing differential corrections to a set of 
initial elements. The results are reported in Table and the 
(O-C) residuals of the individual velocities from this solution 
are listed in Table ^ 

3.3 The adopted spectroscopic orbits 

The agreement between the Tel Aviv and CfA analyses of 
the spectroscopic orbits is satisfactory for all of the orbital 
parameters except the velocity amplitudes for the outer or- 
bit, Ka and Kb, and the center-of-mass velocity, 7, where 



the differences exceed three times the combined internal er- 
ror estimates. These differences lead to a relatively large 
difference in the mass ratio qab = Ka '■ Kb, which we use 
to divide up the total mass from the astrometric orbit in the 
next section. Apparently there is a systematic error lurking 
in one or both of our analyses, but we have not been able to 
identify the culprit. Therefore we have chosen to adopt final 
orbital parameters that are simply the average of the val- 
ues from the two solutions, as reported in Table ^ For the 
final errors we adopted half the difference between the two 
solutions, or the average of the two internal error estimates, 
whichever was larger. 

The agreement between the astrometric orbit and our 
adopted spectroscopic orbit is well within the errors, which 
are much smaller for the astrometric solution. By conven- 
tion, the longitude of periastron for an astrometric orbit 
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refers to the secondary, but for a spectroscopic orbit it refers 
to the primary. After this difference of 180° is taken into ac- 
count, the values for lj agree well. 



4 THE MASSES AND ORBITAL 
INCLINATIONS OF GLIESE 644 

4.1 The masses of the three stars 

The combination of spectroscopic orbits for the inner and 
outer binaries together with an astrometric orbit for the 
outer visual binary enables us to derive masses for all three 
components of Gliese 644. The total mass of the system can 
be derived from the parallax together with the dynamical 
quantity {a^ / P"^) provided by the astrometric solution. The 
individual masses can then be deduced using the two mass 
ratios derived from the spectroscopic solutions. 

For the parallax we note that the Fourth Edition of 
the General Catalog of Trigonometric Parallaxes (van Al- 
tena, Lee & Hoffleit 1995) gives the weighted average of 
16 determinations as TTtng = 0'.'1548 ± 0'.'0006. The Hip- 
parcos catalog (ESA 1997) lists a considerably larger value 
(tthip = 0'.'1742±0'.'0039), which is almost certainly affected 
by the 1.7-yr orbital motion of the visual pair. Soderhjelm 
(1999) took this orbital motion into account in his reanaly- 
sis of the Hipparcos data and derived vr = '.' 1556 ± '.'0018, 
in excellent agreement with the ground-based value. These 
determinations are plotted in Figure |^ including the 16 in- 
dividual ground-based values. A valuable check is provided 
by Gliese 643, which is a common-proper-motion companion 
to Gliese 644 and is therefore expected to be at the same 
distance (see Section 7). The value of its parallax as mea- 
sured by Hipparcos is tthip = 0'.'1540 ± 0'.'0040, again in 
good agreement with the ground-based value for Gliese 644. 
We adopt the ground-based value for the parallax of Gliese 
644 because it has the smallest estimated error. 

Combining (a^/P^) from our adopted astrometric orbit 
with the parallax we get that the total mass is 



Mtot = 



(aVP' 



= 1.049 ± 0.020 (1.9%) Mq . 



(1) 



An alternative way to derive the semi-major axis of the 
wide binary, and consequently the total mass of the sys- 
tem, is to use the radial-velocity amplitudes derived from 
the spectroscopy together with the inclination angle derived 
from the astrometry. This yields Mtot = 1.17 ± O.4OM0. 
Note that the error for this approach is an order of magni- 
tude larger than the one corresponding to the mass derived 
from the astrometry alone, and the total masses from the 
two approaches agree within this large error. 

The mass ratios given by our adopted spectroscopic 
orbits are qb = Meb : Afea = 0.906 ± 0.015 for the in- 
ner orbit and qab = Mb '■ Ma = 1.56 ± 0.17 for the 
outer orbit. Combining these mass ratios with the defini- 
tions Mtot = Ma + Mb and Mb ~ A/pa + Meb, we conclude 
that 



Ma = 



1 



1 + <?AB 



-A^tot = 0.410 ± 0.028 (6.9%) Mq , 



(2) 



gAB 



{1 + qB){l + Qab) 



Mtot = 0.336±0.016 (4.7%) Mq ,(3) 
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Figure 7. Parallax determinations for Gliese 644. The individ- 
ual ground-based results are plotted as circles. The weighted av- 
erage is indicated by the dotted line. The arrow points to the 
corresponding error bar. Soderhjelm's (1999) redetermination for 
Gliese 644 is plotted as a filled square, replacing the Hipparcos 
value plotted as an open square. 



Mb 



qBQAB 



(1 + gB)(l -HgAB 



-Mt, 



0.304±0.014 (4.7%) M0.(4) 



These values are similar to the ones found by 
Soderhjelm (1999), who derived Ma = 0.41 ± 0.04 and 
A^B ~ O.66±O.O6M0, and to the ones derived by Segransan 
et al. (2001), who obtained masses of 0.4155 ± 0.0057, 
0.3466 ±0.0047 and 0.3143 ±0.0040, for the masses of A, Ba 
and Bb, respectively. 



4.2 The orbital incHnations 

We can now go back and calculate the inclinations of the two 
orbits relative to our line of sight, using the spectroscopic 
values for M sin'^ i adopted in Table ^ together with the 
individual masses derived in the previous section. We get 

iAB = 17?4±1?4 or 162?6±1?4, is = 16?3±0?3 or 163?7±0?3 .(5) 

The ambiguity comes from the fact that the radial-velocity 
data can not specify whether the orbital motions are ret- 
rograde or direct. The inclination of the outer orbit from 
the astrometric solution, iAB = 163? 1 ± 1?6, is close to the 
retrograde value from the spectroscopic orbit. 

The inclinations of the two orbits do not provide enough 
information to derive the relative angle between the two 
orbits, 0, which can be written as (Mazeh & Shaham 1976): 

cosff) — coszb cosiAB + siu IB sin ZAB cos(f2B — S^ab) , (6) 

where f^B and Qab are the position angles of the line of 
nodes of the inner and outer orbits, respectively. Since f^B 
is unknown, we can only limit (ji (Batten 1973): 



lA 



JAB < 4> < iA + IAB, 



(7) 



which in our case results in 
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1?1 ± 1?4 < < 33?7± 1?4. (8) 

Equation (8) implies that the relative angle between the 
two orbits could be as large as 34° . However, to make the two 
inclinations come out so close requires a special orientation. 
It is more natural to suppose that the two orbits really are 
almost coplanar. 

5 PHOTOMETRY 

The Hipparcos Catalogue lists V = 9.02 mag for the to- 
tal combined light of Gliese 644, while SIMBAD cites 6 de- 
terminations that lead to V = 9.02 ± 0.01 mag. To divide 
the total light between the three components we use the 
adopted values of the spectroscopic light ratios at 5187 A 
derived by TODCOR, a = 0.569 ± 0.034 and 0.358 ± 0.024. 
After the application of small corrections to move these 
light ratios to 5550 A, the value we adopted for the ef- 
fective wavelength of the V band, these ratios become 
0.576 ± 0.034 and 0.364 ± 0.024. Dividing the V brightness 
with these light ratios, and applying our adopted parallax of 
n = 0'.'1548 ± 0'.'0006, we find that My {A) = 10.69 ± 0.02, 
Mv(Ba) = 11.29 ± 0.05, and Afv(Bb) = 11.79 ± 0.05 mag. 

As a check on the spectroscopic light ratios derived with 
TODCOR, we calculate the total light that those ratios pre- 
dict for the B subsystem and compare it with the magnitude 
difference between A and B found by others. The simple av- 
erage of the 70 visual magnitude differences reported in the 
Washington Double Star Catalog (Worley & Douglass 1996) 
gives 0.103 ± 0.014. Soderhjelm (1999) reports a difference 
in the Hipparcos magnitudes of 0.11 based on his reanalysis 
of the Hipparcos data. These magnitude differences are both 
consistent with our result that the total light of B is fainter 
than A by 0.07 ± 0.05 mag in the V band. 

The three new speckle observations reported in this pa- 
per provide infrared magnitude differences, and show that B 
is brighter tha.n A at J, H, and A" by 0.51±0.01, 0.56±0.02, 
and 0.66 ± 0.06 mag, respectively. The fact that B is fainter 
than A in the visual but brighter in the infrared is qualita- 
tively consistent with the fact that the two components of 
B are cooler than A. 

However, when we use our derived masses to predict 
the magnitudes expected for A and B in detail, a problem 
emerges. This is illustrated in Figure where we plot the 
observed absolute magnitudes for A (tesselated stars) and 
B (filled triangles) at V, J, H, and K, using the appar- 
ent magnitudes for Gliese 644 (all three stars together) of 
J = 5.37, H = 4.67, and K = 4.38. We used Leggett (1992) 
observations and applied small corrections to convert them 
to the Johnson system. We also plot the predicted abso- 
lute magnitudes as a function of wavelength from 5 Gyr 
solar-metallicity models (solid lines, Baraffe et al. 1998) and 
[Fe/H] = —0.5 metal-poor models (dash-dotted lines, Baraffe 
et al. 1997). In both cases the lower line of each pair corre- 
sponds to the prediction for A, and the upper line for B. The 
infrared absolute magnitudes predicted for A are slightly 
brighter than observed, while the corresponding predictions 
for B are typically too faint. In the visual band, the solar- 
metallicity model is slightly too bright for both stars. 

These discrepancies become more obvious when we plot 
the magnitude differences (B— A), as shown in panel b. The 
5 Gyr solar-metallicity and [Fo/H] = —0.5 metal-poor model 




V J H K 

I , , , , I , , , , I I , , : I , , I I 

0.5 1.0 1.5 2.0 2.5 

Wavelenf^Lh (/im) 



Figure 8. Absolute magnitudes (panel a) and magnitude dif- 
ferences (panel b) for Gliese 644 A and B. In panel a the ob- 
served A magnitudes are denoted by stars, and those of B by 
triangles. The predicted absolute magnitudes for A and B from 
5 Gyr solar-metallicity models are shown as solid lines, and from 
[Fe/H] = —0.5 metal-poor models as dash-dotted lines. In panel b 
a 0.1 Gyr solar-metallicity prediction is plotted as a dashed line. 
The plus signs present the magnitude differences predicted by the 
Henry & McCarthy (1993) empirical mass-luminosity relations. 

predictions are again plotted as solid and dash-dotted lines, 
respectively, while a 0.1 Gyr solar-metallicity prediction is 
plotted as a dashed line. Although metallicity and age both 
have some effect on the predicted magnitude differences, the 
effect is much smaller than the discrepancy with the ob- 
served magnitude differences. To show that the models ap- 
pear to be reliable, we plot (as plus signs) the magnitude 
differences predicted by the Henry & McCarthy (1993) em- 
pirical mass-luminosity relations. We are puzzled by these 
discrepancies and see no easy explanation. 



6 THE MASS-LUMINOSITY RELATION AND 
GLIESE 644 

The individual masses that we have derived for the three M 
dwarfs in the Gliese 644 system have formal errors of 5 to 

7 %, which puts these results into a category with only a 
dozen or so stars that have masses smaller than 0.6 Mq and 
mass uncertainties similar or better than ours. At this level 
of accuracy we can begin to make meaningful tests of the 
theoretical models for low-mass stars. 

It has been traditional to assess the mass- luminosity 
relation using the absolute V magnitude versus mass dia- 
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Figure 9. The mass-luminosity relation for stars on the lower 
main sequence with mass determinations more accurate than 5%. 
The three components of Gliese 644 are plotted as tesselated stars. 
For comparison we also show the determinations for Gliese 644 
by Segransan et al. (2001). The components of the three eclipsing 
binaries YY Gem, GJ 2069A, and CM Dra are plotted as filled 
circles, and the components of the astrometric systems Gliese 860, 
Gliese 866, and Gliese 791.2 as open circles. The solid line is the 
solar-metallicity 5 Gyr isochrone, and the dashed line is the metal- 
poor 5 Gyr isochrone with [Fe/H] = —0.5. The corresponding 10 
Gyr isochrones are plotted as dotted lines. 



gram. We can plot all three stars in the Gliese 644 system 
on this diagram, because TODCOR has provided us with 
the spectroscopic light ratios needed to divide up the hght 
between the two components of the unresolved inner binary 
at a wavelength close to the V band. Thus we use the AIv 
versus mass diagram to assess the present status of the con- 
frontation between the observations and the BarafTe et al. 
(1997, 1998) theoretical isochrones for low-mass stars. 

In Table 6 we list the masses and absolute visual magni- 
tudes obtained here for Gliese 644, as well as data for other 
stars below 0.6 Mq with uncertainties in the masses smaller 
than 5%, which we take as a benchmark for comparison. 
These measurements are shown in Figure ^ We did not list 
and plot Gliese 570 BC (Forveille et al. 1999) because that 
close system has not been resolved in the visual, and there- 
fore observed V magnitudes are not available for the two 
components. 

In the figure we plot solar-metallicity 5 Gyr isochrone 
from Baraffe et al. (1998), and metal-poor 5 Gyr isochrone 
with [Fe/H] = -0.5 from Baraffe et al. (1997). The corre- 
sponding 10 Gyr isochrones are also plotted. As can be seen 
from the figure, some of the observed points are far off from 
both theoretical lines. 

Forveille et al. (1999) suggest that the Baraffe et al. 
models (1997, 1998) may be missing some opacity sources, 
and that the solar-metallicity isochrones should be lower 
than they appear in Figure pi This might produce better 
agreement with the observed results for CM Dra, Gliese 866 
and 860, in which case the good agreement of our three 
points for Gliese 644 might imply that the system is slightly 
metal poor and/or very young. Can the much more serious 
discrepancies for GJ 2069A, as suggested by Delfosse et al. 



(1999a), and for Gliese 791.2, be explained by extreme metal 
richness? The answer to this question will have to wait for 
the calculation of theoretical isochrones for super-metal-rich 
models and for the observational determination of accurate 
metallicities for these stars. In the meantime we must con- 
clude that the mass-luminosity relation near the bottom of 
the main sequence is still not well understood, at least in the 
visual band. Note that the above dicussion did not take into 
account the possiblity of stellar spots, which could affect the 
stellar luminosity. In the future, it might be better to move 
to infrared magnitudes, where the effects of stellar spots on 
the photometry are less severe. 

Figure ^ shows that changing the stellar age from 5 to 
10 Gyr has no significant effect on the theoretical isochrones 
for stars less massive than about 0.5 Mq. However, changing 
the metallicity from [Fe/H] = 0.0 to —0.5 has a large effect. 
To allow a meaningful confrontation between the models and 
the observations in the My versus mass diagram, it is nec- 
essary to have accurate metallicities for the stars involved, 
good to 0.1 in [Fe/H] or better. Otherwise, the metallicity 
uncertainties are likely to dominate the uncertainties from 
the photometry, parallaxes, and masses. Unfortunately, ac- 
curate metallicity determinations are still beyond the state 
of the art for stars near the bottom of the main sequence. 
Of all the systems plotted in Figure for example, only 
for CM Dra has the metallicity been analyzed in some de- 
tail, with ambiguous results (Viti et al. 1998). They find that 
CM Dra may be somewhat metal poor, which would be con- 
sistent with its unusually high space motion, or it may have 
almost solar metallicity, but it is unlikely to be metal rich. 

Actually, it is possible that metallicity alone is not 
enough to understand the mass-luminosity relation. For ex- 
ample, it may turn out that it is not valid to assume that 
the helium abundance scales with the metallicity, in which 
case it will also be necessary to determine the bulk helium 
abundance. Stellar age is certainly another factor that can 
change the stellar luminosity. We expect therefore that stud- 
ies of stars in clusters will (continue to) be important for 
progress in testing stellar models, because additional infor- 
mation about the age and metallicity can be derived for 
cluster stars. 



7 THE DISTANT COMPANIONS OF GLIESE 
644 

In this section we discuss two additional stars associated 
with the Gliese 644 triple system. One of them is Gliese 643 
(=Wolf 629, a = 16:55:25.26, S = -8:19:21.3 [J2000], 
V — 11.74 mag), at a projected separation of about 70" 
from Gliese 644. The other companion, at a separation of 
220", is the faint star vB 8 (=Gliese 644C, a = 16:55:35.74, 
5 = -08:23:36.0 [J2000], V = 16.80 mag). The common 
proper motion of Gliese 644 and 643 (Wolf 1919) and their 
similar parallax (e.g., ESA 1997) strongly indicate that 
they are indeed physically connected. Van Biesbroeck (1961) 
found that vB 8 also shares with Gliese 644 the same proper 
motion, attesting to its physical association with the system. 

vB 8 (spectral type M7.0 V) and vB 10 (M8.0 V), a 
companion to Gliese 752 also found by van Biesbroeck, were 
long considered to represent the bottom of the stellar main 
sequence. Although many cooler and less luminous objects 
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Table 6. Lower main-sequence stars with mass and V absolute magnitude better than 5%. 



Star 


Mv 








Mass 






Ref. 


Gls 791.2 A 


13.38 


± 


0, 


.03 


0.2866 




0.0061 


Benedict et al. 2000 


Gls 791.2 B 


16.65 


± 


0, 


.10 


0.1258 


± 


0.0029 




Gls 860 A 


11.81 


± 


0, 


,07 


0.257 




0.011 


Henry & McCarthy 1993 


Gls 860 B 


13.39 


± 


0, 


.06 


0.172 


± 


0.008 


Henry et al. 1999 


GJ 2069 A 


11.7 


± 


0, 


,2 


0.4329 




0.0018 


Delfosse et al. 1999a 


GJ 2069 B 


12.45 


± 


0, 


.2 


0.3975 


± 


0.0015 




Gls 866 A 


15.34 




0, 


.14 


0.1216 




0.0029 


Delfosse et al. 1999b 


Gls 866 B 


15.58 


± 


0, 


,07 


0.1161 


± 


0.0029 




Gls 866 C 


17.34 


± 


0, 


,45 


0.0957 


± 


0.0023 




YY Gem A 


8.99 


-I- 





,08 


0.588 


± 


0.022 


Andersen 1991 


YY Gem B 


8.99 


± 


0, 


,08 


0.588 


± 


0.022 




CM Dra A 


12.75 


± 


0, 


,04 


0.2307 


± 


0.0010 


Metcalfe et al. 1996 


CM Dra B 


12.90 


± 


0, 


,04 


0.2136 


± 


0.0010 




Gls 644 A 


10.69 


± 


0, 


,02 


0.410 ± 0.028 


This paper 


Gls 644 Ba 


11.29 


± 


0, 


,05 


0.336 ± 0.016 




Gls 644 Bb 


11.79 


± 


0, 


,05 


0.304 ± 0.014 





are now known, vB 8 is still of interest, because its mass 
must be very near the substellar limit of 0.08 Mq. Using the 
updated mass- luminosity relation of Henry et al. (1999), we 
derive from the My values of Gliese 643 and vB 8 (12.69 and 
17.75 mag) mass estimates of 0.19 and 0.08 Mq , respectively. 

Joy (1947, see Abt 1973) obtained a few spectra of 
Gliese 643 and reported that its radial velocity was vari- 
able. Since that time Gliese 643 has usually been reported 
in the literature as a spectroscopic binary (e.g., Eggen 1978, 
Johnson 1987). To follow its orbital motion, we secured 83 
spectra of Gliese 643, spread over 5843 days, but could not 
find any significant radial- velocity variation. Our velocities, 
derived using Gliese 725B as the template, have an r.m.s 
of 0.63 kms~^, slightly larger than our typical errors for M 
dwarfs, and the probability is small, but we can not find 
any periodicity in the velocities or any orbital solution that 
pass our usual tests for being significant. If Gliese 643 is 
a spectroscopic binary, the orbital amplitude must be less 
than about 0.5 kms~^. Our individual velocities for Gliese 
643 are reported in Table ^. 

The mean radial velocity we find for Gliese 643 is 
15.81±0.07 kms~^on the CfA system, i.e. with the same zero 
point as the CfA velocities for Gliese 644. This is very close 
to the center-of-mass velocity of Gliese 644, 7 = 15.10±0.21. 
The projected separation between Gliese 644 and 643 is 
about 450 AU, which corresponds to an orbital velocity of 
about 1.5 kms~^ for a circular orbit. The observed velocity 
difference is only 0.7 kms~^, consistent with the interpre- 
tation that Gliese 643 is gravitationally bound in an orbit 
with Gliese 644. 

We have also secured one low S/N spectrum of vB 8 
with the Digital Speedometer on the MMT and found its 
radial velocity to be 15.68 ± 1.09 kms~^, consistent with 
the interpretation that vB 8 is bound to the Gliese 643/644 
system in a hierarchical quintuple configuration. In this pic- 



Table 7. Radial velocities and internal error estimates (km s ^) 
for Gliese 643 (first 20 lines) 



HJD Vr (Tint 



2445780, 


,9993 


16, 


,44 


±0, 


,41 


2446511, 


,9753 


15, 


,84 


±0, 


,26 


2446512, 


,8511 


15, 


,82 


±0, 


,45 


2446513, 


,8675 


15, 


,87 


±0, 


,59 


2446520, 


,0150 


15, 


,94 


±0, 


,33 


2446520, 


,8916 


16, 


,00 


±0, 


,85 


2446523, 


,8416 


15, 


,78 


±0, 


,39 


2446537, 


,9360 


15, 


,73 


±0, 


,20 


2446538, 


,9771 


16, 


,22 


±0, 


,26 


2446539, 


,8486 


15, 


,79 


±0, 


,60 


2446540, 


,8417 


15, 


,03 


±0, 


,32 


2446540, 


,9364 


15, 


,83 


±0, 


,21 


2446541, 


,8916 


15, 


,17 


±0, 


,22 


2446565, 


,7301 


15, 


,42 


±0, 


,43 


2446568, 


,7910 


15, 


,97 


±0, 


,34 


2446569, 


,6356 


15, 


,75 


±0, 


,79 


2446569, 


,7865 


15, 


,58 


±0, 


,53 


2446569, 


,7964 


16, 


,27 


±0, 


,59 


2446608, 


,6348 


15, 


,54 


±0, 


,94 


2446612, 


,7619 


14, 


,80 


±0, 


,69 



ture vB 8 moves around the center of mass of the Gliese 
644/643 system. 

The projected separation between vB 8 and Gliese 644 
is only three times larger than the projected separation be- 
tween Gliese 643 and Gliese 644. Such a small ratio usually 
renders triple systems dynamically unstable (e.g., Eggelton 
& Kiseleva 1995, Kiseleva et al. 1995). Therefore, we suggest 
that the actual separation between vB 8 and Gliese 644 is 
significantly larger than its projected separation, by at least 
a factor of two. 
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8 ON THE FORMATION AND EVOLUTION 
OF THE SYSTEM 

Our picture of the Gliese 644/643/vB 8 system can be sum- 
marized as follows: 

• All five stars reside in an hierarchical system. The or- 
bital sizes are of the order of 2000 (projected), 500 (pro- 
jected), 1 and 0.05 AU for vB 8, Gliese 643, Gliese 644AB, 
and Gliese 644B, respectively. 

• The masses of vB 8, Gliese 643, Gliese 644A/Ba/Bb are 
of the order 0.1, 0.2, 0.4, 0.3, 0.3 Mq, respectively. 

• The innermost two orbits, in the Gliese 644 system it- 
self, are probably coplanar. 

In this section we briefiy discuss two extremely sim- 
plistic scenarios for the formation and evolution of the sys- 
tem. First we consider a scenario in which the quintuple 
was formed by a sequence of fragmentation events during 
the coUapse of a molecular cloud core (e.g., Burkert & Bo- 
denheimer 1993) leading directly to the hierarchical con- 
figuration now observed. The cloud core collapsed until it 
was a few thousand AU in diameter, at which point it frag- 
mented into three parts. The two smaller parts ended up 
as vB 8 and Gliese 643. The larger fragment continued to 
contract by another three orders of magnitude, until a hy- 
drostatic core was formed with a radius of a few AU (e.g., 
Larson 1969). When the temperature of the core reached 
about 2000 K the molecular hydrogen began to dissociate, 
resulting in a second collapse within the hydrostatic core 
(e.g. Larson 1969). During the second collapse the core may 
have been able to fragment again into the components of 
the visual binary, Gliese 644A and Gliese 644B, the latter 
of which fragmented again into Ba and Bb. 

This scenario is not free of some problems. Although 
fragmentation calculations have successfully simulated the 
formation of wide binaries during the collapse of molecular 
cloud cores, there is still some uncertainty about whether 
or not this step can succeed in forming a stable binary with 
separation as small as 1 AU (e.g.. Boss 1989, Bonnell & Bate 
1994, Bate 1998). Moreover, there have not yet been success- 
ful fragmentation calculations for the formation of very close 
binaries, such as Gliese 644B with a separation of only 0.05 
AU. Whatever mechanism was responsible for the formation 
of Gliese 644B system, it is interesting that the mass ratio 
for this binary is close to unity, a mass ratio which is ap- 
parently quite frequent in spectroscopic binaries (Tokovinin 
2000). A mechanism that can lead, under the right condi- 
tions, to nearly equal masses in a close binary has been stud- 
ied by Bate & BonneU (1997) and Bate (2000). Initially, two 
stellar seeds form, with masses at most a few percent of the 
final stellar masses. Most of the mass of the collapsing cloud 
core is still in a gaseous envelope surrounding the seeds. The 
accretion of material onto the seeds then builds the masses 
up to the final stellar values. The smaller of the two seeds 
is favored in the amount of mass that it receives, and this 
tends to equalize the masses of the final components of the 
binary, after the envelope material has been exhausted. 

If the two orbits in the Gliese 644 system are coplanar, 
as suggested by our spectroscopic results, then this might 
lend support to fragmentation as the formation mechanism 
for the innermost binary; it is natural to suppose that the 
two orbital planes for Gliese 644 should both be oriented per- 



pendicularly to the angular momentum vector of the molec- 
ular cloud core from which they formed. 

Dynamical evolution in a small-N cluster (e.g., McDon- 
ald & Clarke 1993, 1995) is another possible scenario for 
the formation of the Gliese 644/643/vB 8 system. In this 
scenario all five stars were formed independently within 
a loosely self-gravitating small-N cluster of stars. Grav- 
itational interactions between the cluster members then 
evolved the orbits of the five stars, the ones we see now, into 
a hierarchical configuration, while ejecting the other stars 
from the cluster. In this scenario, the more massive stars are 
expected to settle into more tightly bound orbits, while the 
less massive stars are raised into less tightly bound orbits, 
just as is observed for the system. However, this scenario 
does not easily account for the coplanarity of the Gliese 644 
triple, a feature often found in other triple systems as well 
(Fekel 1981, Tokovinin 1997). 

Another possible way to explain the coplanarity of the 
close triple Gliese 644 involves tidal interactions. It is well 
known that a large relative inclination between the two or- 
bits of a triple can induce strong oscillations of the inner 
eccentricity (e.g., Mazeh, Krymolowski & Rosenfeld 1997), 
which can make the system dynamically unstable. There- 
fore, the low relative inclination of the two orbits in the 
Gliese 644 triple might be the result of some selective evo- 
lutionary process. Triple systems with large relative inclina- 
tions simply could not survive. This argument still has to 
be worked out, though, because the amplitude of the eccen- 
tricity modulation does not depend linearly on the relative 
inclination (e.g., Holman, Touma & Tremaine 1997), and 
there are medium-sized angles, on the order of 10-20 de- 
grees, that do not produce large eccentricity modulations. 

A potentially interesting feature of the Gliese 644 triple 
— the very marginal evidence for a minute variation of the 
outer orbit's eccentricity — might also be a result of tidal 
interactions (Mazeh & Shaham 1979, Paper II). A few spec- 
troscopic triple systems have aheady been observed to dis- 
play direct or indirect evidence for long-term tidal modula- 
tions (Mayor & Mazeh 1987, Mazeh 1990, Jha et al. 1997), 
of which HD 109486 is the best example (Jha et al. 2000). 
Right now the effect observed in Gliese 644 is much too 
marginal to rely on. Many more years of observations will 
be needed to see if the effect is real. 

The Gliese 644/643/vB 8 system is the nearest known 
quintuple stellar system. The next nearest system with five 
or more objects is Castor (=Gliese 278) at 16 pc (Tokovinin 
1999). In order to better understand how multiple systems 
of low-mass stars form and evolve, it would be extremely in- 
teresting to find more M-star multiple systems in the solar 
neighborhood and to compare their features. In this vein, it 
is amazing that another triple M-star system exhibits very 
similar features to Gliese 644. This is Gliese 866, which was 
recently studied by Delfosse et al. (1999b) and by Woitas et 
al. (2000). Delfosse et al. derived the two orbital periods, of 
803 and 3.78 days, and masses of 0.12, 0.11 and 0.096 Mq. 
Woitas et al. (2000) derived a period of 821 days for the 
outer orbit. Although the masses in the Gliese 644 triple are 
almost three times larger, in both systems there are three 
M stars with similar masses, where the most massive star 
is 30% more massive than the lightest one, and the period 
ratio is slightly larger than 200. Moreover, both systems are 
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probably very close to coplanar. One wonders if these fea- 
tures are common in M-star triples. 
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